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ABSTRACT 

The  dynamics  of  the  sidewall-type  surface  effect  ship 
(SES)  for  vertical  plane  heave  motion  in  response  to  waves 
are  established  through  theoretical  model  and  experimental 
investigation.       The  theoretical  models  are  digitally  simu- 
lated and  compared  to  an  experimental  tow  tank  model  SES 
(length  to  beam  ratio  of  6.5)  for  verification.       The  exam- 
ination of  transient  and  frequency  responses  provides 
insight  into  the  fundamental  dynamic  characteristics. 
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NOMENCLATURE 

A  Waterplane  area  of  cushion,    I  »b. 

A  Exit  air  flow  area. 
e 

A   ,  ,  A  Bow  and  stern  exit  air  flow  areas, 

eb'     es 

A  Entrance  air  flow  area. 
o 

A  Sidewall  waterplane  area,    I  •  y 
sw  ^  »         Jsw 

a  Added  mass. 

B  Beam  overall. 

b  Cushion  width. 

b  Cushion  damping  coefficient. 

b  Sidewall  damping  coefficient. 

bn  b     +  b 
1  c         w 

C   ,  C  Quadratic  orifice  coefficients. 
e      o 

F.  External  applied. force. 

g  Gravity  constant. 

h  Exit  gap  height. 

h  Equivalent  cushion  height  at  equilibrium. 

h.  Equilibrium  gap  height. 

k  Polytropic  constant. 

k1  Sidewall  stiffness  coefficient. 

k2  Equivalent  cushion  stiffness  coefficient. 


k  .      Equilibrium  buoyancy  force  on  immersed  sidewall,    2pV 
L  Length  of  SES  overall. 

I  Length  of  cushion. 

:M  Air  mass  in  cushion  region. 


* 


M,  Supply  flow  rate. 


'°A11  equilibrium  conditions  are  indicated  by  a  zero  subscript. 


NOMENCLATURE  (continued) 

M.  Air  mass  flow  rate  into  cushion, 
in 

M  Air  mass  flow  rate  out  of  cushion, 
out 

m  Mass  of  the  craft. 

p  Ambient  air  pressure,    absolute. 

3. 

"~r>  Cushion  air  pressure,    absolute. 

*p  Supply  air  pressure  from  fan,    absolute. 

5 

t  Time. 

u  Craft  forward  velocity. 

*V  Cushion  volume. 
c 

V,  Duct  volume, 
d 

V  Volume  of  immersed  sidewall  at  equilibrium, 
sw  n 

co  Wave  frequency. 

co  Encounter  frequency. 

CO  Natural  frequency. 

x  Horizontal  position  along  craft  length. 

y  Sidewall  thickness. 

z  Wave  height  at  any  point. 

z  Vertical  height  of  craft. 

z  Cushion  height. 
c  & 

z   .  Water  height  inside  cushion  at  bow. 

z  Water  height  inside  cushion  at  stern, 
os  6 

z  Sidewall  height. 

z  Average  wave  height  inside  cushion. 

£  Damping  ratio. 

rl  Wave  height. 

'"All  equilibrium  conditions  are  indicated  by  a  zero  subscript. 


a 
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NOMENCLATURE  (continued) 


X  Wave  length. 

p  Water  density. 

p  Air  density. 


p  Cushion  air  density. 

p  Supply  air  density. 
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CHAPTER   I 
INTRODUCTION 

1 .  1       Background 

In  the  past  few  years  substantial  interest  has  developed  in  surface 
effect  ships  (SES),    specialized  waterborne  air  cushion  vehicles  which 
can  provide  high  speed  waterborne  transportation.       The  SES  configura- 
tion being  studied  is  of  the  captured  air  bubble  type.       It  employs  the 
air  cushion  vehicle  principles,  utilizing  a  pressurized  compartment  of 
trapped  air  to  support  the  majority  of  the  craft's  weight.      A  smaller 
fraction  of  the  total  support  force  comes  from  the  hydrostatic,   hydro- 
dynamic  and  aerodynamic  forces.       The  pressurized  compartment  of  air 
(cushion)  is  contained  by  rigid  sidewalls  along  the  craft  length  and  by 
seals  both  fore  and  aft.       See  Fig.    1.1. 

In  the  ocean  environment,   the  principal  vertical  plane  motions  of 
interest  from  the  aspect  of  seakeeping  and  habitability  are  heave,    surge, 
and  pitch.       So  as  to  provide  reasonable  estimates  of  these  motions, 
effort  is  directed  towards  the  development  of  tools  that  will  provide 
information  on  the  expected  behavior  of  the  full-scale  craft.       The 
principal  tools  used  here  are  experimental  model  tests  and  mathematical 
simulation  via  digital  computer. 

The  dynamic  simulations  are  helpful  where  known  Froude  scaling 
laws  do  not  hold.       Since  this  is  true  of  the  heave-pressure  coupled 
motion  for  the  SES,    then  the  simulation  model  becomes  a  very  useful  tool 
for  full-scale  motion  prediction.       This  does  not  preclude  the  use  of 
experimental  scale  models  but,    rather,    it  means  that  once  digital 
simulations  accurately  predict  experimental  model  results,   then  this  new 
tool  can  be  easily  adjusted  for  known  parameter  changes.       Therefore,    it 
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13 
can  predict  more  accurately  full-scale  craft  dynamics. 

1.  2       Scope  of  Thesis 

The  principal  thrust  of  this  thesis  initially  leads  in  two  directions. 
These  are  (1)  to  develop  a  mathematical  model  of  the  equations  of  motion 
for  digital  simulation  and  (2)  to  build  and  test  a  scale  model  of  an  SES  in 
the  MIT  tow  tank  facility.       Attention  is  focused  to  the  single  degree  of 
freedom  in  heave  only.       Heave  motion  of  the  SES,    being  the  main  contrib- 
utor to  the  craft  vertical  acceleration  levels,    is  a  good  starting  point  for 
digital  simulation  and  experimental  verification. 

Chapter  II  explains  the  modeling  of  the  fluid  suspension  system. 
The  primary  components  to  the  model  are  developed  for  later  use  in  the 
simulation  Models  A,    B,    and  C  found  in  Chapters  III  and  IV  where: 

Model  A  is  the  SES  craft  transient  performance  analysis  on 
cushion  with  no  entrance  or  exit  air  flow.       (Principally  this  is  an 
inverted  open-bottom  container  suspended  by  a  pressurized  air  mass  on 
the  water  surface.  ) 

Model  B  is  the  same  as  Model  A  except,    now,    air  flow  is  considered. 

Model  C  is  the  simulated  responses  to  regular  waves. 

The  most  basic  model  consists  primarily  of  the  following  major 
components:     (1)  a  cushion  region  formed  under  the  craft  that  is  main- 
tained by  low  pressure  air  for  craft  support;     (2)  a  nonflexible 
sidewall  structure;     (3)  the  forward  and  after  seals;    and  (4)  a  fluid 
source  for  feeding  air  to  the  cushion. 

The  simulated  transients  and  wave  frequency  analysis  are  to  be 
compared  to  experimental  model  tests  for  verification.       Comparisons  of 
heave  amplitude,    phase  and  acceleration  levels  are  illustrated,    along 
with  their  time  histories  at  selected  frequencies.       The  experimental 


14 


model  tests  were  conducted  at  zero  forward  speed. 

1.  3       Prior  Work  Related  to  Heave  Motions 

Historically  there  has  been  interest  in  the  analysis  of  air  suspen- 
sion systems  in  heave  from  two  major  technologies:    (a)  the  land-based 
air  cushion  vehicles,    and  (b)  the  sea-based  air  cushion  craft.       It  is  the 
latter  field  to  which  this  thesis  addresses  itself,    although  useful  informa- 
tion and  analysis  has  been  obtained  from  both  areas. 

The  most  recent  and  extensive  analysis  and  simulation  studies  for 
SES  craft  have  been  conducted  by  Oceanics  [1].'      Simulation  model 
analysis  of  the  SES  vertical  plane  dynamics  has  been  performed,   though 
verification  is  not  extensive  in  the  report. 

Other  land-based  air  cushion  suspension  studies  have  been 
conducted  for  the  Department  of  Transportation  at  MIT  [2  and  3],    the 
results  of  which  have  been  widely  published. 


[  ]  indicate  references. 


Note:      From  henceforward,    "model"  and  "modeling"  shall  refer  to  the 
theoretical  model  for  simulation.       In  order  to  avoid  confusion  in 
terminology,    "experimental  model"  or  "SES  experimental  craft" 
is  used  when  referring  to  the  tow  tank  test  craft. 
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CHAPTER    II 
THE  HEAVE  SUSPENSION  MODEL 

2.  1       The  General  Suspension  Configuration 

A  nonlinear  dynamic  model  for  heave  motion  of  a  rigid  body  is 
developed  for  the  configuration  illustrated  in  Fig.    1.1.       The  basic  sus- 
pension system  consists  of  (1)  cushion  region,    (2)  supply  system, 
(3)  seal  system,    and  (4)  sidewalls.       In  the  configuration  a  fan  supplies 
flow  through  a  constant  volume  duct  and  orifice  restriction  to  the  cushion 
chamber  (plenum).       The  flow  from  the  plenum  escapes  through  the 
sealing  gap  area  formed  between  the  seal  and  the  water  surface.       The 
fluid  suspension  is  basically  a  three  energy  port  device  where  energy 
transfer  between  the  waves,    air  suspension  and  craft  takes  place  as 
illustrated  in  Fig.    2.  1.       In  developing  models  for  these  parts  of  the 

Gravity 
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Craft 


Forces 

and 
Motions 


Suspension 
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Water  Waves 
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-^    Acceleration 
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Fig.   2.  1      Craft  -  Suspension  -  Wave  System. 
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system,   the  following  assumptions  are  made: 

(1)  The  cushion  and  source  absolute  pressures  are  comparable 
to  the  absolute  ambient  pressure.       For  the  experimental 
tow  tank  model  the  cushion  pressure  is  on  the  order  of  one- 
tenth  psig  and  the  supply  pressure  of  one  to  two  psig. 

(2)  The  external  flow  and  pressure  fields  induced  around  the 
vehicle  while  it  is  moving  with  a  forward  velocity  are  not 
considered  to  affect  the  exit  flow  from  the  cushion.       This 
assumption  is  not  valid  at  high  speeds,    especially  at  the 
bow  seal.       Compensation  for  this  effect  could  be  modeled 
into  the  system  for  more  extensive  study. 

(3)  The  cushion  chamber  base  area  A     is  constant. 

(4)  The  pressure  p     and  density  p    in  the  cushion  region  are 
spatially  uniform. 

(5)  The  Froude-Krylov  hypothesis  assumes  the  wave  system 
is  not  appreciably  altered  by  the  presence  of  the  craft. 
Thus,    it  is  assumed  that  the  wave  is  not  diffracted  or 
attenuated  as  it  passes  under  the  cushion  chamber. 

(6)  The  craft  is  assumed  to  be  a  rigid  body. 

The  principal  forces  acting  are  identified  in  the  rigid  body 
equations  of  motion  for  heave,    Fig.    2.  2. 

dz- 

(m+a)  -5-i-     =    Z+Z+Z+Z  (1) 

■     '    dt  c         w  s  g 


where 


Z       =    forces  due  to  the  cushion  region 
c  b 

-     "Ab  tec-Pa*  "Vl 


(la) 
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(a)  Coordinate  System. 

Craft  Weight 


Hydrostatic  and  Dynamic  Wave  Forces 
(b)  Cushion  Forces. 

Craft  Weight 


Hydrostatic  and  Dynamic 
Sidewall  Forces 


(c)  Sidewall  and  Craft  Forces 


Stern  Seal 


Hydrostatic  and 

Dynamic  Forces 
Bow  Seal 


(d)  Bow  and  Stern  Seal  Forces. 


Fig.    2.2      Coordinates  and  Principal  Heave  Forces. 
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where 


Z         =      sidewall  forces 
w 


•k, (z. -z    )  -  b    z1    -  k, 
1     1      w  w    1  1 


-      mg+F1 

a      =      added  mass 
A,       =      waterplane  area  of  cushion  base 

A  =     waterplane  area  of  the  sidewall 

sw  r 

b  =  cushion  damping  coefficient 

b  =  sidewall  damping  coefficient 

F  =  external  force  acting  on  the  SES 

k  =  stiffness  coefficient  of  sidewall 

=      2p  A 

K      sw 

k1       =      vertical  buoyant  force  of  immersed 

sidewall      =      2  p  V 

sw 

m  =  mass  of  craft 

p  =  ambient  pressure  absolute 

p  =  cushion  pressure  absolute 

z1  -  craft  position 

z         =      average  external  water  elevation 

w  & 

V  =      volume  of  immersed  sidewall  at 

sw 

equilibrium 


(lb) 


Z         =      forces  due  to  seal  (considered  small  ,,    . 

s  (lc) 

compared  to  other  external  forces) 
Z         =      gravity  and  externally  applied  force 

o 


(Id) 


19 


Therefore,    substituting  (1,  a-d)  into  (1 ), 


(m+a) 


dz 


■A 


b(Pc"Pa)  "  \z1  -  k^-zj  -  kj  +  mg  +  F1 


(2) 


where 


b    +  b 
c        w 


In  the  analysis  and  later  in  the  experiment,  the  forces  acting  on 
the  seals  were  assumed  to  be  relatively  small  compared  to  the  cushion 
forces. 


mg  +  F. 


External 

Water 

Surface 


Internal 

Water 

Surface 


Cushion 
Region 


2.  2       The  Cushion  Region  Model 

Within  the  plenum  chamber,    conservation  of  mass  is  required; 
thus,   the  inlet  supply  flow  (M.    )  less  the  flow  through  the  exit  gap  (M      .) 
must  equal  the  rate  of  change  of  mass  within  the  chamber: 

dM 


dt 


M.     -  M      . 
in  out 


(3) 


M 


air  cushion  mass 


The  pressure  in  the  plenum  chamber  may  be  related  to  the  mass 
(M   )  and  chamber  volume  (V   )  by  the  equation  of  state.       A  polytropic 
relationship  is  assumed  for  the  cushion  chamber. 


(  )  indicate  equations. 
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p        =      (M    /V   )k  constant  (4) 

The  polytropic  process  is  felt  to  be  between  isothermal  (k  =  1.0)  and 
isentropic  (k  =  1.4). 

The  total  cushion  volume  V     at  any  instant  in  time  may  be 
expressed  as  the  product  of  the  base  area  A,    and  the  sum  of  the  plenum 
equilibrium  height  h   ,    craft  elevation  z.,    and  the  water  surface  elevation 


V 


where 


V        =      Afe(-z1  +  hQ  +  zQ)  (5) 


h        =      equivalent  cushion  height  at  equilibrium 
z         =      average  wave  height  under  cushion 
z        =      static  depression  of  water  plus  the  external 
wave  elevation  =  (p   -p   )/p    +  z 


2.  3       Cushion  Flow  Entrance  and  Exit 

•  • 

From  (3)  the  entrance  flow  M.     and  exit  flow  M      .  are  modeled  as 

in  out 

follows.       First,   for  purposes  of  the  experimental  craft  verification  the 

inlet  mass  flow  M.     is  considered  constant  since  the  entrance  flow  came 

in 

from  a  high  pressure  source  where  the  dynamics  of  the  cushion  region 
did  not  affect  the  mass  flow  to  any  degree.       Secondly,    in  modeling  an 
actual  system  where  fan  and  duct  dynamics  are  important,   the  mass  flow 
into  the  cushion  is  modeled  according  to  the  classical  quadratic  orifice 
equation  for  non-compressible  flow.       Non-compressible  flow  may  be 
assumed  for  the  absolute  source  and  cushion  pressures  since  they  are 
comparable  to  the  absolute  ambient  pressure  (p   -p     =  p    -p      ^  1  psi). 

S         cL  C         c.l 

Mass  flow  into  the  cushion  is.- 
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M,-„      -      CA(2(po"PjPs^  (6) 


in  O      O         '  S     *C     "Si 


The  mass  flow  that  exits  from  the  plenum  must  either  reverse  and 

pass  through  the  entrance  orifice  in  the  case  where  p    >p     or  the  flow 

passes  under  the  bow  and  aft  seals.      Again  the  M        is  modeled 

to  out 

similarly  to  the  entrance  mass  flow  with  the  exception  that  now  the  exit 
leakage  area  A     is  varying. 


e 


Mout      =      CeV^c-Pa^s/g*'  (?) 


where 


A        =      entrance  flow  area  to  cushion 
o 

A        =      variable  exit  flow  area  under 
e 

seals 


C    ,  C  =      discharge  coefficients 

(assume  .   0.  6) 
p        =      fan  supply  pressure 


2.  4       Exit  Air  Flow  Region 

The  exit  area  is  the  opening  between  the  seal  and  the  water  surface 

through  which  the  relatively  low  pressure  air  escapes  at  the  bow  and 

stern.       This  assumes  that  the  sidewalls  allow  no  leakage.       The  total 

leakage  area  is 

A        =     A  ,    +  A  (8) 

e  eb         eg 


where 


A        =      total  exit  area 
e 


A  ,  ,    A  =      bow  and  stern  leakage  areas 

eb'       es  s 

The  bow  and  stern  exit  areas  are  related  geometrically  to  the  craft 
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seals  and  water  vertical  position  by: 


eb 


es 


b(-z,   +  z   .    +  h.  )        if  >0 
1         ob        1 


if  £0 


b(-z,   +  z       +  h,  )        if  >0 
1         os         1' 


if  £0 


(9) 


(10) 


where  the  wave  height  at  the  bow  and  stern  respectively  are  given  by: 


'ob 


os 


cushion  water  height  at  bow 
cushion  water  height  at  stern 
equilibrium  gap  height 


2.  5       The  Fluid  Source  and  Feed  System 

This  section  on  the  fluid  source  is  presented  for  future  detailed 
studies  where  the  duct  and  fan  dynamics  may  influence  the  full-scale 
craft  results.       The  fan  and  duct  dynamics,    as  given  here,    were  not 
needed  for  the  verification  of  the  experimental  model  tests.       The  use  of 
a  constant  mass  flow  rate  source  in  the  experimental  test  apparatus  made 
it  feasible  to  focus  attention  directly  to  the  fluid  suspension  (cushion) 
system.       The  feeding  system  which  supplies  air  flow  to  the  cushion  in 
the  full-scale  craft  consists  of  a  fan  and  a  feed  duct.       The  fan  may  be, 
in  general,    assumed  to  operate  at  a  constant  supply  pressure.       This 
assumption  is  valid  if  the  source  pressure  remains  nearly  equal  to  the 
fan  operating  point  equilibrium  pressure.       However,    for  large  deviations 
from  the  operating  point,    the  change  in  supply  pressure  as  a  function  of 


supply  flow  should  be  considered. 

The  fan  curves  identified  in  Fig.    2.  3  may  be  represented  by  the 
nonlinear  functional  relation- 
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Mf.  =    f(Ps-Pa) 


(11) 


Operation  about  the  equilibrium  point  may  be  linearized  for  simulation 
and  analysis  wherein  the  mass  flow  rate  IVL  is  given  as.- 


M. 


M 


fo 


so 


+  2 


(12) 


where 


M.      =      supply  flow  rate 
equilibrium  supp 
equilibrium  supply  pressure 


M~      -      equilibrium  supply  flow  rate 


rso 

The  supply  feeding  duct  is  assumed  to  act  as  a  fluid  capacitance 
with  duct  resistance  lumped  into  the  plenum  inlet  supply  area  resistance. 

The  equation  of  state  for  a  polytropic  pressure  density  relation- 
ship similar  to  (4)  combined  with  the  conservation  of  mass  may  be 
written  as; 


VdPs       • 
kp  ps 


M,  -  M. 
f  in 


(13) 


where 


p       =      supply  air  density 


V 


duct  volume 


2.  6       External  Disturbances 

The  transient  performance  study  for  Models  A  and  B  considers  an 
instantaneous  step  increase  in  force  F..   applied  to  the  model  craft.       In 
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the  later  study  of  Model  C  the  external  disturbance  comes  from  the 

waves  z     .       The  external  wave  disturbance  acts  to  change  the  air  volume 
w  b 

and  hydrostatic  force  on  the  sidewalls. 

The  typical  seaway  considered  in  the  analysis  is  a  regular  sinus- 
oidal wave  of  height  TJ     and  frequency  co.       In  the  SES  coordinate  system, 
the  origin  is  considered  at  the  longitudinal  center  of  gravity.       Therefore 
the  regular  wave  elevation  for  head  seas  is  described  by: 


Stern 


-.6/2 


Bow 


+1/2 


X 


where 


V 


o 
\ 

u 


CO. 


CO 
X 


.     ,2lTx  .. 

T?    sm(  — —  +  co  t) 
o  X  e 


length  between  seals 
wave  height 
wavelength  of  wave 

craft  forward  velocity 

2      , 
encounter  frequency  =  cot  +  co     u/g 

for  head  seas 

wave  frequency 

position  along  craft  length 

from  midship  positive  forward 


(14) 
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The  instantaneous  wave  height  along  the  length  of  the  cushion  and 
at  the  seals  is  important  in  the  accurate  model  representation  of  the  real 
system.      In  the  first  case  the  water  height  at  the  bow  and  stern  seals 
acts  as  an  external  disturbance  that  changes  the  equilibrium  gap  height 
h     caused  by  air  flowing  from  the  cushion  under  the  seals.       The  wave 
height  at  the  bow  and  stern  seals,   without  taking  into  account  the 
depression  due  to  cushion  pressure,   is  given  by: 

z   .       =      77    sin(to  t  +  TT^/X)  (15) 

ob  'o  e  ' 

z  =      77    sin(o)  t  -  TTi/X)  (16) 

os  'o  e  ' 

The  average  wave  height  under  the  cushion  is  determined  by 
integration  along  the  cushion  length.       The  average  wave  height  under  the 
cushion  is  given  by: 


z 

w 


Therefore 


w  $_  77 


77      ri/2 

■—  \  sin(2rrx/X  +  co  t)dx 

S.    J.^/2  e 


7?   X 

£—  sin(7Tj2/X)  sin  to   t  (17) 


This  model  is  representative  as  long  as  the  air  pressure  inside  the 
cushion  is  considered  a  constant  throughout  and  there  is  no  localized 
wave  pumping. 

2.  7       Combined  Model  Equations 

The  system  equations  defining  the  heave  dynamics  of  the  craft  for 
Models  A,    B,    and  C  in  Chapters  III  and  IV  are  given  below.       This 
includes  the  equations  from  sections  2.1   -  2.5. 

(2)    dz    /dt  =  Zj 

(m+a)dz1/dt  =  -b^  -k%1(z1-zw)  -Ab(pc-pa)  -kj  +  mg  +  Fj 
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(3)    dJVI   /dt    =     M.     -  Mf 


(4)    p      =    (M   /V   )k  constant 


(5)    V       -    A^(-  z,   +  h     +  z    ) 
c  b         1         o        o 

where  z     =  (p     -  p   )/p  +  z 


(6)  M.       =    C   A   (2(p     -  p   )p  /   )2 

in  o    o      *s       *c  's/g 

(7)  M      .     =    C    A   (2(p     -  p   )p    I  ) 

out  e    e      ^c      ^a  ^c/g 


(8)    A       =    A  .    +A 

e  eb        es 


(9)    A  .     =        b(-z.   +  z   ,   +  h. )         if  >0 
eb  r  1         ob         1  - 


(•  1  OD  1     -| 

I  0  J 


if  so 


(10)    Aes    =        b(-z1+zos+hl)        if>0 


(-  1  US  I     ^ 

i  0  J 

(11)    Mf    =    f(pg  -  pc) 

(13)   (vdps/kpso)Ps         u. 


if  £0 


(18) 


in 

2.  8       Linearized  Model  Equations 

A  linearized  set  of  equations  for  the  heave  mode  was  developed  and 
used  along  with  the  nonlinear  equations  (18)  for  determining  effects  of 
cushion  parameters.       The  linearized  equations  are  given  in  state 
variable  form  for  Models  A,    B,    and  C.       Excluded  from  these  three 
models  are  the  duct  and  fan  dynamics  since  the  experimental  SES  was 
supplied  from  a  constant  mass  flow  source. 


d 


— ■               — 

A  z^ 

A  z 

= 

A  M 
c 

, 

AAA 
21       A22      ^23 


A31      ° 


A  Z- 

23 

A  z. 

+ 

33 

A  M 
c 

,.,_ 

— 

0  0 


B21     B22 


0  B 
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A  F. 


A  z 


w 


(19) 
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Ap, 


Cll      ° 


13 


A  z- 

A  z 

+ 

A  M 
c 

0  D 


12 


A  F 


A  z 


w 


where 
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k33 


B 


21 


322 


B 
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'11 


'13 


D 


12 


-(k^  +  k2)/(m+a) 

-b1/(m+a) 

-AbC13/(m+a) 

-K3Cn  -  K4(-l  +Cn/p) 

-K3C13  -  K4C13/p 

l/(m+a) 

(kj   -  AbD12)/(m+a) 

-K3D12-K4(1  +Dia/p) 

&PCA>Z!  =  P/(h0P/kPCo+  !) 

bp    /5M       =    (h    p/M      )/(h    p/kp       +1) 
*c'  C  O^'       CO'       o^'     rco 


bp    /bz 
^c'       w 


-(p/(h0p/kpco  +  l) 


where 


K3      '      Co^c//  Ae/[2(Pco  "  Pa^ 
K4      "      Co(2Pc/g>*    2b<Pc-Pa^ 


AbCll 
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CHAPTER  3 
SURFACE  EFFECT  SHIP  TRANSIENT  RESPONSE 

3.  1       Introduction 

The  basic  equations,    as  developed  in  Chapter  2  and  summarized  in 

2.  7  and  2.  8,    are  considered  in  the  study  of  the  transient  response  to 
external  inputs.       The  theoretical  results  of  various  responses  are  com- 
pared to  experimental  data  for  verification.       Two  theoretical  models 
that  closely  repeat  the  experimental  tests  are  designated  by  the  titles, 
Model  A  (without  air  flow),    and  Model  B  (with  air  flow).       As  a  result  of 
the  transient  responses,    the  natural  frequency  and  damping  of  the  SES 
are  obtained. 

Prior  to  reading  in  detail  the  analysis  in  Chapters  3  and  4,    turn  to 
Appendix  A  for  discussion  of  the  experimental  model  and  testing. 

3.  2       Transient  Performance  Parameters 

In  the  transient  performance  analysis,    linearized  cushion  model 
performance  is  compared  to  experimental  test  data.       Response 
characteristics  are  determined  for  several  types  of  disturbances  which 
are  explained  in  2.6. 

The  Model  A  equations  are  representative  of  the  heave  dynamics  of 
the  SES  on  the  air  cushion  without  air  flow  passing  through  the  boundaries 
of  the  craft.       Basically,    the  Model  A  dynamic  equations  are  second 
order  wherein  the  craft  vertical  position  z     and  velocity  z     are  the 
primary  state  variables.       The  cushion  pressure  and  other  desired  out- 
puts are  dependent  variables  of  this  system.       The  results  of  this 
analysis  are  the  determination  of  the  craft  natural  frequency  and  damping. 
The  experimental  data  transients  are  conducted  while  the  craft  remains 
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stationary  in  the  horizontal  plane. 

From  (18),    Model  A  nonlinearized  equations  are  given  with  the 
following  conditions; 

dM   /dt,    dM./dt,    dp   /dt,    M.   ,    M     .     -    0 
c'     '         f     '      ^s'  in'       out 

Upon  linearization,    the  resultant  equations  in  state  variable  form 
are  listed  in  (19)  where: 

A23'    A31'    A33'    B32'    C13     =    ° 

The  second  order  system's  characteristic  equation  in  linear  form 
is; 

(m+ajSj   +  bjZj   +  (^  +  Abp/(hop/kpco  +  \))zl    =    0  (20) 

The  quadratic  equation  indicates  a  natural  frequency  for  the  coupled 
heave-pressure  mode  given  by: 

k.   +  A    /(h    /kp       +  1)- 


00      = 

n 


■1  b' 


°         co  1  (21) 

+  a  J 


m  + 
The  damping  ratio  (relative  to  the  critical  damping)  is  given  by 

£    =    b1/(2oon(m+a))  (22) 

In  (20)  the  numerator  terms  are  a  measure  of  the  system's  stiff- 
ness.      The  first  term  k..   is  the  stiffness  coefficient  of  the  sidewalls. 
The  second  group  of  terms  of  the  numerator  is  the  stiffness  coefficient 
due  to  the  air  cushion.       It  should  be  noted  that  the  theoretical  stiffness 
coefficient  is  approximately  equal  to  the  waterplane  area  of  the  cushion 
and  sidewalls  times  the  water  density. 

All  parameters  of  the  second  order  system  can  be  determined 
directly  (see  Table  A.  1)  except  for  the  damping  b1    and  added  mass  a. 
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Estimates  of  these  two  variables  may  be  derived  from  comparison  of  the 
transient  heave  response  of  the  model  and  experiment.       As  a  check  on 
these  results  from  theory  [5],    estimates  of  the  damping  and  added  mass 
for  oscillating  pressure  fields  on  a  free  surface,    as  formulated  by 
Stoker  (1957),   were  used.       Curves  for  added  mass  and  damping  from 
[5]  are  found  in  Appendix  B.       Note  that,    in  order  for  correlation  between 
experiment  and  theory  to  be  valid,   the  seal  and  sidewall  effects  must 
also  be  estimated  since  19%  of  the  total' waterplane  area  is  due  to  the 
sidewall. 

Figure  3.  1  is  the  response  for  a  step  input  in  force  acting  on  the 
craft.       This  force  was  37.  5  lb,    equal  to  100%  of  the  craft  weight.       The 
natural  frequency  and  critical  damping  ratios,    calculated  by  comparing 
the  digital  response  and  experimental  data,    are  approximately: 

to       =9.2  rad/sec 

n  ' 

4     =    0.3 

where  m  =  1.16,    a  =  4.34  (lb  sec2  /ft),   kj  +  kg  =  461  lb/ft,  bj  =30  lb  sec/ft. 

Determination  of  the  added  mass,    a,    from  theory  based  upon  a 
two-dimensional  oscillating  pressure  field  [5]  was  in  agreement  with  the 
experimental  results.       Theory  from  [5]  (see  Appendix  B),    indicates 
that  the  damping  ratio  due  to  the  cushion  is  on  the  order  of  £  =  0.  3. 
This  appears  to  agree  well  enough  with  the  experimental  results,   though 
some  discrepancy  due  to  the  sidewalls  is  certain  to  be  present. 
Reference  [l  ]  indicated,    though,    that  the  damping  due  to  the  pressure 
field  is  negligible,   which  is  contrary  to  what  [5]  and  experimental  results 
indicated. 

In  this  study  it  was  expected  that  there  would  be  more  of  a 
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difference  between  theory,    based  solely  upon  the  cushion  dynamics  [5], 
and  the  experimental  results.      This  is  because  the  experimental  craft 
had  solid  sidewalls  and  seals  which  covered  19%  of  the  total  waterplane 
area  affected.       But,   within  the  accuracy  of  the  experimental  results,    no 
detectable  difference  was  found.       Some  theoretical  estimates  of  the 
hydrodynamic  force  terms  for  the  sidewalls  could  be  found  through  the 
use  of  slender  body  theory,    but  for  this  thesis  no  such  estimates  are 
made.       Future  studies  should  look  into  this. 

Typical  velocity,    acceleration  and  pressure  transients  corres- 
ponding to  Model  A  are  shown  in  Fig.    3.2. 

3.  3  Model  B  Transient  Response 

The  transient  response  analysis  for  Model  B  is  performed  in  the 
same  manner  as  in  Model  A.       The  principal  addition  to  the  dynamic 
equations  is  the  entrance  and  exit  mass  flows  described  in  (3).      The 
entrance  mass  flow  is  assumed  to  be  a  constant  while  the  exit  mass  flow 
is  described  by  the  quadratic  orifice  equation  (7)  where  the  exit  area  and 
cushion  pressure  vary  with  time. 

Again  the  resultant  natural  frequency  co     and  damping  ratio  £  are 

determined.       The  state  variables  of  the  third  order  system  are  heave  z., 

velocity  z        and  air  cushion  mass  M   .       The  resultant  output  variables 

that  are  observed  are  the  cushion  pressure  p     and  acceleration  z    .      The 

1  r c  1 

model  equations  from  (18)  have  the  following  conditions  imposed; 

dM   /dt,    dp    /dt    =    0      and      M.       =    const, 
c'  s'  in 

Upon  linearization,   the  resultant  equations  in  state  variable  form 
are  listed  in  (19).       These  equations  are  used  in  the  transient  response 
study. 
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Prior  to  comparing  model  simulations  to  test  data,    some  observa- 
tions of  the  experimental  tests  are  made. 

3.  4  Transient  Response  from  Experimental  Tests 

The  SES  experimental  model  response  with  air  flow  through  the 
seals  presented  some  interesting  observations.       Figure  3.  3  shows  two 
transient  response  curves  for  two  test  conditions  nearly  identical  except 
for  the  seal  design.       Since  the  seals  used  in  the  experiments  were  non- 
flexible,   the  mass  air  flow  rate  would  be  expected  to  vary  more  dramat- 
ically here  than  for  a  design  using  a  flexible  seal  that  contours  the  water 
surface. 

In  test  1  of  Fig.    3.  3,   the  heave  response  shows  a  flattening  at  the 
top  of  the  curve  when  the  craft  rises  to  0.  02'  above  equilibrium.       In 
test  2  the  craft  rises  to  0.038'  above  equilibrium.      Why  the  difference? 


In  test  1  the  seal  design  was  such  that  a  vertical  gate  seal  of 


i  n 


thick  plexiglass  was  used  to  contain  the  cushion  at  the  bow  and  aft.  In 
this  case  the  exit  gap  length  was  very  short,  resulting  in  a  more  rapid 
bleed-off  of  air.       In  test  2  the  seal  was  much  thicker  (2  inches).       This 
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Seal  design 
for  test  1. 


Seal  design 
for  test  2. 
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seal  is  faired  to  allow  the  water  to  flow  smoothly  around  it  without 
creating  excess  turbulence.       The  air  escaping  under  the  seal  has  to 
travel  a  larger  distance  before  escaping  from  the  plenum.       For  this 
reason  the  air  loss  is  less,    and  the  mean  pressure  level  inside  the 
cushion  is  higher. 

As  a  result  of  the  seal  design,    the  heave  response  of  test  2  has  a 
higher  overshoot  than  in  test  1.       The  corresponding  pressure  drops  in 
tests  1  and  2  were  8.6  and  7.2  lb/ft",    respectively.       Therefore,   the 
seal  design  can  be  expected  to  influence  the  response.       All  further 
experimental  tests  use  the  thick  seals. 

Another  important  observation,  which  is  not  in  the  simulation,  is 
the  standing  wave  that  is  observed  being  created  by  the  seals.  This  has 
the  effect  of  causing  local  fluctuations  in  cushion  volume  and  the  exit  air 
gap  other  than  the  normal  rise  and  fall  of  the  water  level. 

3.  5  Comparison  of  Simulation  and  Experimental  Transients 

The  transient  response  curve  for  the  linearized  equations  is  com- 
pared to  the  test  data  in  Fig.    3.  4.       The  natural  frequency  for  the  linear- 
ized equation  of  co     =9.2  rad/sec  is  the  same  as  in  Model  A.       The 
^  n 

curves  of  lower  damping  are  representative  of  the  real  situation  where 
the  faired  seal  (test  2)  is  used.       This  has  a  damping  ratio  on  the  order 
of  4  =  0.  3.       The  experimental  response  (test  1)  has  a  £  =  0.  4  for  the 
first  quarter  cycle  of  the  nonlinear  response.       The  addition  of  non- 
linearity  in  the  exit  gap  area  causes  more  overshoot,    with  the  same  value 
of  §  =  0.  3,    than  was  found  in  the  totally  linearized  model  equations.       The 
value  of  £  =  0.  3  is  carried  through  to  the  frequency  response  analysis. 


38 


r- 

CO 

i 

CD 

U 

U 

o 

fn 

r- ( 

CO 

rt 

on 

u 

o 

>H 

cu 

o 

CM 


co 


T-H 

G 

cu 

•  iH 

co 

CO 

o 
II 

o 
ii 

o 

G 

cd 

H 

i — i 

jjj»     _ 

T— 1 

cd 
C 

C          C         N 

CM 

0 

CO          CO            « 

Eg03 

0i 

••-I 

■S               -rH                  M 

ii 

u 

f- 

f- 

1                  _ 

cu 

0)           CD           -C 

c      a       a 

3C 

a 

CO 

X 

5 

i* 

1 — 1 

i — i 

W 

£ 

£ 

O 

CD 
T3 
O 

o 

T3 
G 

=8=         =tfc         rf 

^_ 

cd 

-t- 
Cf 

a 

est 
ear 

cd 

cu 

o 

cu 
en 

pq 

. — i 

cu 

E- 

E- 

G 

•rH 

CD 

0) 

-a 
o 

i 

1 

J 

J       " 

o 

a 

S 

i-i 

^-* 

•rH 

1 
J 

T 
i 
i 
i 

H 

cd 

0) 
G 

^ 

• 

o 

CM 

• 

o 

o 

g 

o 

CO 

•rH 

Sh 
cd 

a 

S 
o 
U 

• 
CO 

tub 

•r-t 
ft 

o 


o 
I 


CO 

o 


o 

i 


CM 
O 


o 
I 


o 
I 


o 
o 


o 
o 


CM 

o 


CO 

o 


o 


o 


*"       c 

>       o 
cu 

n 


39 


Another  result  of  the  Model  B  is  that  the  addition  of  the  air  mass  flow 
has  little  effect  upon  the  natural  frequency  and  critical  damping  when 
compared  to  Model  A  without  air  flow.       This  is  a  very  likely  result 
since  the  time  constant  for  filling  the  air  cushion  region  is  long  com- 
pared to  the  wave  frequency. 

The  last  comparison  to  be  made  is  the  acceleration  and  pressure 
responses,    Fig.    3.  5,    of  the  digital  simulation  and  experimental  tests. 
Agreement  appears  to  be  quite  good  except  for  the  higher  frequency 
oscillations  which  are  thought  to  come  from  the  bubbling  action  of  the  air 
through  the  water  around  the  seals.       Also  some  impact  forces  may  be 
caused  by  the  water  slapping  the  seals.       These  localized  dynamic 
responses  are  not  accounted  for  in  the  Model  B. 
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CHAPTER  4 
SES  FREQUENCY  RESPONSE  TO  REGULAR  WAVES 

4.  1       Introduction 

With  the  conclusions  of  the  transient  response  the  following  chapter 
presents  the  frequency  response  analysis.       Heave  motion  is  the  principal 
dynamic  response  to  be  analysed,    since  it  can  be  easily  compared  to 
experimental  data.       Information  on  acceleration  and  cushion  pressure  is 
also  presented.       The  behavior  in  regular  waves  was  measured  in  the  tow 
tank  with  the  SES  craft  constrained  to  heave  only  and  operating  at  zero 
forward  velocity.       The  simulation  model  is  presented  first  and  then  com- 
pared to  the  experimental  results. 

4.  2       Model  C  Frequency  Response 

The  modeling  of  the  heave  dynamic  frequency  response  is  primarily 
based  on  linear  theory  with  certain  nonlinear  effects  also  included.       The 
linearized  equations  of  Model  C,    as  developed  in  Chapter  2,    are  included 
in  (19).       The  third  order  system  equations  are  the  same  as  given  in 
Model  B  except  for  the  following  additions: 

1.  The  external  disturbances  considered  are  the  wave  pumping  (17) 
and  the  wave  elevation  at  the  bow  and  stern  seals,    (15)  and  (16). 

2.  The  principal  nonlinearity  that  is  included  in  the  systems 
equations  is  the  gap  exit  area  (9),    (10)  and  (11),    wherein  the  leakage  area 
is  allowed  to  vary  as  a  function  of  the  craft  seal  elevation  and  wave 
height.       The  model  considers  a  nonflexible  seal  attached  to  the  bow  and 
stern.       With  the  craft  stationary  in  the  x-y  plane,    air  leakage  exits  from 
both  seals.       When  the  craft  is  moving  forward  then  the  principal  leakage 
is  from  the  stern  seal  only. 
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3.     The  body  suspension  forces  acting  on  the  craft,    aside  from  the 
pressure  force,    represent  the  hydrostatic  and  hydrodynamic  effects  due 
to  body-fluid  interaction  (including  wave-induced  effects)  for  the  sidewalls. 
Estimates  for  added  mass  and  damping  come  from  the  transient  analysis 
and  [5]. 

4.  3      Heave  Suspension  Performance 

In  order  to  illustrate  the  nature  of  the  craft  responses  to  waves, 
digital  computer  simulations  and  experimental  time  histories  are 
presented.       The  results  of  typical  simulations  appear  in  Fig.    4.1.       The 
responses  to  regular  waves  illustrate  the  magnitude  and  phase  of  the 
heave,    acceleration  and  pressure  to  the  input.       The  constraint  on  the 
leakage  area  given  by  A,    ^  0  is  presented  to  illustrate  the  basic  nature  of 
the  outputs.       All  motions  are  quite  regular  and  the  variation  in  gap  area 
is  the  only  nonlinearity.       The  behavior  of  the  vertical  CG  acceleration  is 
approximately  proportional  to  but  slightly  lagging  the  pressure  as  the 
wavelength  increases.       The  heave  mean  position,   though  not  shown  in 
Fig.    4.  1,    shifts  downward  as  the  gap  area  increases  until  a  new  equilibrium 
position  is  attained. 

Figure  4.  2  illustrates  typical  time  histories  of  the  motion  found 
during  the  experiments  in  the  tow  tank.       Several  comparisons  between 
the  simulations.    Fig.    4.1,    and  experiments,    Fig.    4.2,    are  made: 

1.  Note  the  nonlinearity  in  the  experimental  heave  response.       The 
air  leakage  at  the  bow  and  stern  appears  to  be  causing  the  majority 
of  the  nonlinearity.       The  simulation  model  of  the  air  gap  leakage 
is  more  linear. 

2.  The  experimental  heave  motion  shows  significant  flattening  of 
the  response  near  the  peak,    which  does  not  appear  in  the  linear 
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Fig.    4.  1      Simulated  Heave  Time  Histories  of  SES  to  Regular  Waves. 
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Fig.    4.2      Experimental  Heave  Motion  Time  Histories  for  SES. 
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simulations. 

3.     The  heave  motion  lags  the  wave  input  and  is  fairly  consistent 

between  experiment  and  simulation. 

The  pressure  and  acceleration  time  histories  from  several  experi- 
mental results  are  presented  in  Fig.    4.  3.       The  high  frequency  oscilla- 
tions in  the  pressure  and  accelerations  are  felt  to  be  caused  by  the  air 
escaping  under  the  seal  due  to  the  bubbling  action  which  transmits  pres- 
sure pulses.       Likewise,    they  appear  in  the  acceleration.       The  high 
frequency  oscillations  occur  when  the  craft  is  reaching  its  peak  in  heave 
and  the  air  exit  gap  opens.       These  high  frequency  oscillations,   being  due 
to  local  effects  around  the  seals,    should  be  given  careful  attention  in 
future  studies.       The  high  frequency  oscillations,    though,    appear  to  have 
little  effect  in  the  overall  heave  motion  prediction  when  comparisons 
between  simulation  and  experiment  are  made  on  the  macro  level.       In 
addition,    the  solid  seals  attached  to  the  craft  are  transmitting  the  forces 
due  to  the  wave  impact  on  the  seals.       A  flexible  seal  should  result  in 
lower  wave  impact  force. 

In  order  to  give  a  better  picture  of  what  is  taking  place  overall,    let 
us  analyse  the  heave  dynamics  over  a  range  of  wavelengths.       The  heave 
motion  (rms),    Fig.    4.  4,    from  the  digital  simulation  indicates  a  behavior 
somewhat  similar  to  that  indicated  by  the  limit  value  given  in  (2  3),    with 
near  zero  response  at  values  of 

"-  K I  -  # ^  (¥)      .  <23> 

XI H  -  1.  0,   -j,  -i-  and -j.       Therefore,   the  minimum  heave  response  will 

occur  when  the  wave  pumping  is  also  a  minimum,    and  likewise  the  maxi- 

9      2  2 

mum  response  occurs  close  to  kjl    =    2,  ^,   -g-,    and  j.       The  small  differ- 
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ence  between  the  simulation  and  the  values  in  (23)  is  due  to  the  effect  of 
the  natural  heave  frequency  which  is  near  Xj i  =  0.  38. 

The  large  response  characteristics  at  short  wavelengths,    which 
are  due  to  the  direct  influence  of  wave  pumping  as  well  as  to  the  reson- 
ance effect,    result  in  relatively  large  heave  acceleration  (rms)  values  as 
indicated  in  simulation  results,    Fig.    4.  5.       This  condition,    which  is 
generally  known  to  occur  in  air  cushion  vehicles,    is  often  referred  to  as 
"cobblestone"  oscillation.       It  is  evident  when  considering  the  operation 
of  these  craft  over  short  wavelengths.       This  effect,    though,    was  not 
verified  from  the  experimental  data,    but  is  recommended  for  future 
verification. 

The  results  of  the  heave  motion  and  acceleration  from  simulations, 
Figs.    4.  4  and  4.  5,    agree  with  the  simulation  studies  in  [1]. 

With  the  basic  results  of  the  heave  motion  presented,    the  simula- 
tions are  now  compared  to  the  experimental  test  results.       The  heave 
motion  frequency  response  data  from  several  series  of  experimental 
tests  are  summarized  in  Fig.    4.  6.       Data  for  the  response  to  waves  of 
length  A,  =  1.  5  -  12  ft,    and  height  r\     -  2  in,    indicate  that  the  model 
exhibits  the  same  general  response  as  the  experimental  results.       This 
is  •         the  minimums  and  maximums  along  with  the  increasing  heave 
response  amplitude  with  wavelength  are  similar.       Experimental  data 
for  wavelengths  3.1   -  4.6  were  not  available,    Fig.    4.6;    therefore,    the 
results  are  not  conclusive  in  this  region. 

The  most  interesting  result,    Fig.    4.  6,    is  the  shift  in  the  experi- 
mental heave  response  minimums  as  compared  to  the  model.       The  two 
primary  minimums  are  at  A/ 1  =  0.  35  and  1.  47.       At  the  lower  wave- 
lengths another  minimum  appears  near  A/ 1  -  0.22. 
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In  an  effort  to  account  for  the  differences  between  experimental 
and  model  results,    several  ideas  were  put  forward  based  upon  observa- 
tions of  the  cushion  waves.       With  the  lack  of  data  for  X  =  3.  1  -  4.  6,    the 
first  assumption  concerning  the  response  in  this  area  was  to  assume  the 
heave  response  followed  a  smooth  curve.       This  would  then  indicate  that 
the  primary  minimum  heave  response  shifts  to  the  left  and  right  of  the 
model  values  X  =  0.  5  and  1.  0. 

In  order  to  determine  what  may  cause  this  shift,    observations  of 
the  wave  profile  photographically  indicated  some  differences  between  the 
external  and  internal  waves.       Further  experimental  testing  was  per- 
formed by  mounting  wave  probes  internally  and  externally  to  the  SES 
craft.     Experimental  results  for  various  wavelengths,    Fig.    4.  7,    indicate 
that  the  incoming  wave  is  being  altered  in  several  ways.       First  the 
internal  wave  is  lagging  the  external  wave  by  20-30  degrees  at  the  lower 
wavelengths.       Secondly,    secondary  waves  are  appearing.       Thirdly,    the 
internal  wave  is  amplified  at  data  values  A  =  3.  12',    6.  5',    and  9',    and 
slightly  attenuated  at  the  other  wavelengths.       This  indicates  that  the 
seals  and,   to  a  lesser  degree,    the  sidewalls  are  creating  an  internal 
standing  wave.       For  this  reason,    the  minimum  heave  response  would  be 
expected  to  shift.       The  amount  of  shift  would  depend  upon  whether  the 
standing  wave  or  the  external  wave  is  dominant.       The  standing  wave, 
having  its  own  characteristic  eigenvalue,    could  cause  a  shift  in  the 
system's  natural  frequency  response  away  from  the  original  model  values. 

There  is  some  speculation  that  a  dip  in  the  heave  response  might 
occur  in  the  region  A  =  3.1  -  4.  6  where  experimental  data  were  not  avail- 
able,   Fig.    4.  6.       If  this  were  true,   then  the  results  would  indicate  that 
the  minimums  have  shifted  to  the  right  of  model  values.       Again,   the 
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Cushion  Waves 
Ew      =      External  Waves 
Abscissa  -  Time,    1  sec  =  25  div. 
Ordinate    -  1  in.    =  10  div. 

Fig.    4.  7      Comparison  of  Cushion  and  External  Waves. 
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internal  standing  wave  would  be  the  main  contributor  to  such  a  shift. 

Since  no  further  model  simulations  were  performed  that  could 
account  for  the  wavelength  shift  and  experimental  verification  of  the 
heave  response  at  the  missing  wavelengths,    this  is  left  for  future  inves- 
tigation.    Additionally,    investigation  with  the  craft  moving  at  a  forward 
velocity  is  needed.       Sufficient  questions  have  arisen  that  will  hopefully 
induce  others  to  further  the  research. 

Additionally,  data  comparisons  of  heave  response  to  wave  heights 
of  1,  2  and  3  inches,  Fig.  4.  8,  indicate  that  heave  damping  increases 
with  higher  wave  heights  for  A/ 1  <  1.0.  The  maximum  and  minimum 
heave  response  to  waves  of  varying  height  appears  to  shift  toward  the 
model  values  as  the  wave  height  increases,  Fig.  4.  8.  This  might  be 
expected  in  that  the  external  wave  dominates  any  internal  standing  wave. 

This  concludes  the  experimental  tests  conducted  and  the  extent  of 
verification  of  the  simulation. 
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CHAPTER  5 
SUMMARY  AND  RECOMMENDATIONS 

5.  1       Summary 

Theoretical  model  simulations  and  experimental  craft  tests  for  the 
one  degree  of  freedom  heave  motion  of  a  surface  effect  ship  have  been 
conducted.       The  theoretical  third  order  linearized  model  dynamical 
equations  consider  the  following: 

1.  A  constant  air  mass  flow  into  the  cushion. 

2.  Quadratic  exit  gap  orifice  restriction  for  air  flow  out  of 
the  cushion. 

3.  Conservation  of  momentum  of  the  SES  rigid  body  in  the 
vertical  mode. 

4.  Conservation  of  air  mass  inside  the  cushion. 

'5.     Polytropic  relationship  between  pressure,    air  mass  and 
air  volume. 

6.  Stationary  model  in  x-y  plane. 

7.  External  force  inputs  and  regular  wave  disturbances  are 
the  principal  transients  considered  in  analysis. 

In  order  to  simplify  analysis  and  to  verify  the  system's  responses, 
linear  analysis  has  been  developed  and  the  basic  third  order  characteristic 
equation  derived.       A  study  of  the  linear  model  characteristics  under 
external  force  inputs  has  been  conducted  for  the  SES  (X/b  =  6.  5)  with 
constant  air  flow.       In  the  best  fit  of  the  transient  response  simulation  of 
the  model  and  the  experimental  data,    the  natural  frequency  and  critical 
damping  ratio  of  9.2  rad/sec  and  0.3,    respectively,    were  determined. 
From  theory  based  upon  Stoker's  formulation  [5],    the  natural  frequency 
and  damping,    based  upon  an  oscillating  free  surface  pressure  field,    were 
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also  found  to  be  in  good  agreement  with  the  experimental  results.       The 
transient  response  characteristics  of  Model  B  (with  air  flow)  did  not 
appreciably  alter  the  response  of  the  craft  when  compared  to  Model  A 
(without  air  flow).       This  was  due  to  the  time  constant  for  filling  the 
cushion  being  an  order  of  magnitude  greater  than  the  system's  natural 
heave  frequency.       The  experimental  results  verified  this.       The  value  of 
damping  was  determined  to  be  the  critical  measure  in  estimating  the 
amplitude  of  the  transient  response.       As  would  be  expected,    the  increase 
in  damping  helps  to  decrease  the  effects  of  external  disturbances. 

The  outcome,    due  to  variations  in  seal  design,    was  noted  where  the 
effective  exit  gap  length  under  the  seal  influences  the  air  loss.       The 
solid  seal  design  maximizes  the  effects  of  the  nonlinear  characteristics 
of  the  exit  gap  area,    which,    in  turn,    influences  the  heave  motion, 
acceleration  and  cushion  pressure.       These  results  would  likely  be 
attenuated  with  flexible  seals. 

The  sensitivity  of  vehicle  performance  to  wavelength  and  wave 
height  were  studied  with  the  following  primary  results: 

1.  Prediction  of  heave  amplitudes  from  theory  compared  favorably 
with  the  experimental  results.       The  relative  amplitude  to  wave  height 
ratio  was  found  to  decrease  with  increasing  wave  height  for  X/JL  s  1.5  in 
the  experimental  tests.       This  is  indicative  that  damping  increases  as  the 
wave  height  increases. 

2.  The  predictions  of  minimum  wave  pumping  from  theory, 
assuming  the  Froude-Krylov  hypothesis,    should  occur  near  X/Z  =  1,   -j, 
and-q-.       Experimental  evidence,    though,    indicates  that  at  low  speeds  the 
primary  wave  is  altered  by  the  seals  and  sidewalls. 

3.  As  a  result  of  seal  and  sidewall  effects  upon  the  waves  under  the 
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cushion,   the  points  of  minimum  and  maximum  heave  motion  are  altered. 
Experimental  evidence  shows  that  the  minimum  wave  pumping  is  shifted 
to  the  left  and  right  of  the  two  primary  theoretical  values.       For  this 
particular  test  craft,    the  minimums  occurred  at  Xji  =  0.  38  and  1.  47. 
Comparisons  of  the  inside  and  outside  wave  profiles  from  experimental 
tests  indicate  a  general  lag  in  the  wave  under  the  craft  when  compared 
to  the  outside  wave.       Also  the  amplitude  of  the  waves  under  the  cushion 
is  attenuated  in  some  wavelengths  while  others  are  amplified  at  other 
wavelengths.       These  effects  are  not  considered  in  the  simulation  model 
and  therefore  account  for  the  major  difference  between  theory  and  experi- 
ment. 

4.     The  sensitivity  of  acceleration  levels  to  external  disturbances 
was  studied  in  the  simulation  model.       The  accelerations  are  180  degrees 
out  of  phase  to  the  heave  motion.       The  acceleration  (rms)  response  over 
a  range  of  wavelengths  in  the  simulation  gives  the  "cobblestone"  appear- 
ance with  relative  magnitudes  increasing  at  shorter  wavelengths.       No 
experimental  tests  were  conducted  to  verify  this  result. 

Turning  now  to  the  experimental  results,   high  frequency  oscilla- 
tions of  the  cushion  pressure  and  accelerations  were  observed.       This 
phenomenon  occurred  when  the  exit  gap  area  was  opening  and  closing. 
The  air  appears  to  escape  under  the  seals  with  bubbling  action.       This 
effect  is  not  apparent  in  the  simulation  models  since  a  smooth  opening 
and  closing  of  the  exit  area  is  considered.       These  high  frequency 
oscillations,    though,    appear  not  to  alter  the  basic  craft  heave  response. 
Also,   there  may  be  some  water  impact  forces  on  the  solid  seals  which 
could  account  for  some  of  the  high  frequency  accelerations. 

In  conclusion,    the  third  order  linear  model  equations  are  fair 
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representations  of  the  model  dynamics  of  the  craft  and  air  cushion.       The 
principal  differences  in  heave  response  between  theory  and  experiment 
would  appear  to  be  attributable  to  the  water  dynamics  under  the  cushion. 
The  next  step  in  the  modeling  must  allow  for  this  seal-water  interaction 
and  the  resultant  waves  under  the  cushion. 

5. 2       Recommendations 

1.  A  detailed  study  of  the  waves  inside  the  cushion  along  with  their 
cause,   i.  e. ,    seals  and  sidewalls,    and  their  effect  upon  the  cushion 
dynamics  is  needed. 

2.  Modeling  techniques  as  a  result  of  1  could,    therefore,    predict 
the  system  dynamics  to  a  greater  degree  of  accuracy. 

3.  Regions  of  applicability  of  the  Froude-Krylov  hypothesis  should 
be  established. 

4.  Although  the  quadratic  orifice  equations  appear  to  predict  fairly 
accurately  the  air  mass  flow  from  the  cushion,    measurements  of  the  time 
lags  induced  by  the  air  flow  around  the  seals  and  their  effects  upon  the 
overall  dynamics  should  be  considered  in  a  more  detailed  study. 

5.  A  study  of  the  contribution  to  damping  and  added  mass  by  the 
sidewalls  would  be  a  valid  consideration  not  studied  in  this  thesis. 
Slender  body  theory  is  suggested  as  one  approach. 

6.  Further  experimental  model  tests  should  be  conducted  with  the 
craft  moving  with  a  forward  velocity. 

7.  Acceleration  amplitudes  between  theory  and  experiment  need 
verification. 

8.  After  completing  the  heave  study,    the  pitch  motion  of  the  SES 
may  be  investigated. 
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APPENDIX  A 
EXPERIMENTAL  FACILITIES  AND  TESTS 

A.  1      Introduction 

The  SES  with  length  to  beam  ratio  of  6.  5  was  selected  for  study  in 
the  experimental  program  of  this  thesis.       Since  little  open  and  detailed 
literature  is  available  on  the  motion  (heave)  characteristics  of  this 
particular  SES,    a  test  program  was  formulated  and  carried  out.       Results 
from  the  static,   transient,    and  frequency  responses  were  used  in  the 
simulation  study. 

A.  2      Test  Program 

For  the  purpose  of  this  thesis,    the  experimental  test  program  was 
designed  to  accomplish  measurements  of  (1)  the  static  stiffness  coefficient 
(displacement/load  ratio)  of  SES,    (2)  the  transient  response  to  a  step  input 
in  force,    and  (3)  the  frequency  response  of  the  SES  to  a  series  of  regular 
waves. 

A.  3      Test  Craft  and  Facility 

Figure  A.  1  is  a  photograph  of  the  SES  craft  in  the  tow  tank  with  the 
tow  carriage  and  associated  instrumentation.       The  basic  test  craft 
dimensions  and  equilibrium  testing  conditions  are  given  in  Table  A.  1. 
The  test  craft's  sidewall  and  deck  were  constructed  from  ^-inch  plexi- 
glass.      The  sidewalls  and  seals  at  the  bow  and  stern  were  shaped  (faired) 
from  styrofoam  (Fig.    A.  2)  and  attached  to  the  test  craft.       Note:    the 
seals  are  inflexible  and  positioned  for  a  particular  equilibrium  state. 
One  advantage  of  building  the  test  craft  from  plexiglass  is  the  visual 
observation  of  the  free  surface  inside  the  cushion.       Photographs  of  the 
wave  profile  were  taken  from  the  side. 
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Fig.   A.  1      SES  Tow  Tank  Model 
Attached  to  Carriage. 


Fig.   A.  2    SES  Bow  and  Seal  Configuration. 
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TABLE    A.  1 
SES  MODEL  CHARACTERISTICS 

Cushion  Length,    4 (ft) 6.26 

Cushion  Beam,    b(ft)     0.96 

Length  Beam  Ratio 6.5 

Length  Overall,    L(ft)    6.  66 

Beam  Overall,    B(ft)    . 1.19 

Cushion  Area,    A  (ft2)     6.  00 

Sidewall  Height,   z  (ft)     0.  32  8 

Cushion  Height,    z    (ft)     0.  25 

Sidewall  Thickness,   y       (ft)     0.115 

J  sw 

Vertical  Center  of  Gravity,    (ft)    0.386 

Longitudinal  Center  of  Gravity,    (ft)     3.  27    from  stern 

Model  Weight,    (lb) 37.  5 

Specific  Weight  Density,    (lb/ft3)    1.0 

Air  Supply  Flow  Rate    #1  29  ft.  3/min    (STP) 

#2  19ft.3/min    (STP) 

Air  Density 0.  075  lb/ft3 

Water  Density     62.4       lb/ft3 
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Pressurization  of  the  cushion  was  accomplished  by  two  pressurized 
supply  lines,   which  were  3/8-inch  in  diameter,   located  near  the  bow  and 
stern.       Figure  A.  3  illustrates  the  air  supply  and  instrumentation 
involved.       The  air  flow,    hovering  height,    and  craft  weight  may  be 
adjusted  independently.       Nonequilibrium  conditions  were  simulated  by 
introducing  regular  waves.       Several  pressure  transducers,   wave  probes, 
a  heave  rod,    an  accelerometer  and  an  external  flow  meter  were  provided 
to  accomplish  quantitative  studies.      A  chart  recorder  was  available  for 
recording  each  response. 

A.  4      Static  Response 

The  static  response  for  craft  height  versus  load  was  conducted  by 
varying  the  craft  weight  and  measuring  the  vertical  displacement  from 
the  heave  rod.       See  Fig.   A.  4  for  results.       The  observed  stiffness  was 
approximately  linear  over  the  operating  range.       The  linearized  spring 
coefficient  (k     +  k   )  is  equal  to  435  ±10  lb/ft. 

A.  5      Dynamic  Response 

Transient  responses  were  performed  by  removing  weight  from  the 
craft  and  allowing  it  to  reach  a  new  equilibrium.       This  is  equivalent  to  a 
step  input  in  force.      Transient  performance  data  for  heave,    acceleration, 
and  pressure  were  obtained. 

The  frequency  response  tests  were  performed  by  passing  waves  of 
various  lengths  and  heights  past  the  stationary  SES.       Response  data  were 
recorded  for  heave,    acceleration,   pressure,    and  the  wave  profile  inside 
and  outside  of  the  cushion.       Due  to  the  noticeable  difference  between  the 
free  surface  inside  and  outside  the  cushion,   two  different  means  of 
measuring  this  response  were  devised.       Initially,   photographs  were 
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Fig.   A.  4      Experimental  Stiffness  Coefficient  of  SES. 
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taken  of  the  inside  and  outside  waves  through  a  viewing  window  from  the 
side  of  the  tow  tank.       From  these  photos  it  is  possible  to  observe  the 
variance  in  the  waves,   but  quantitative  data  analysis  requires  better 
photography.      The  second  means  of  obtaining  the  time  history  of  the 
cushion  wave  profile  was  by  attaching  wave  probes  at  midship  to  the 
sidewall,   both  inside  and  outside.       The  wave  profile,    relative  to  the 
craft,   was  obtained. 

Discussion  of  the  results  of  these  tests  is  found  in  the  main  text. 
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APPENDIX  B 
ADDED  MASS  AND  DAMPING  FOR  OSCILLATING  FREE  SURFACE 

Figures  B.  1  and  B.  2  were  used  to  estimate  the  cushion  added  mass 
and  damping. 
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APPENDIX  C  -  COMPUTER  PROGRAM  FOR  MODEL  C  68 

c     heave  oymamjcs  or  surface  effect  shif 

C      DIGITAL  SlMyLATICN  USING  DYSYS  (SUBRTUTINF) 

C  vARjApLtS  LISTFO  AP  FOLLOWS 

C  A=  ADoeD  MacS  LB  SEC2/?T 

C  AB«  CUSHION  ARpA  (FT**P) 

C  AE=EUUIL»  AID  rXlT  AREA  FFO"  SEALS  (FT**?) 

C  A  E  B ■  EXfT  ARrA  aT  BOM  (FT**?, 

C  AES=  EXit  AREA  aT  STERN  (  F  T  *  »  2  ) 

C  B«  DAmPIN'G   (|b*SEC/FT) 

C  CO>  CE  B  ENTRANCE  AND  EXIT  cESISTANCE  COEFt   (0.6) 

C  FR=  Alp  TEMP  (r) 

C  G=  GRAVITY  o?,?  ft/SEC*? 

C  He-  EQijILIBPiUm  HEIGHT  OF  CUSHION  (FT) 

C  hi=  CONSTANT  FrR  EQUILIBRIUM  CF  EXIT  A I F  GAP  UN[  Er  SEAL  (FT) 

c  pa*  ambient  aIo  pressure  abs  (lb/ft?) 
c  pc=  Cushion  Pressure  aps(l?/ft2) 

C   HK=  POlyTRORIC  CONSTANT   <1.2) 

C   CJA*  MAj5S  MR  FlOW  «ATE  INTO  CUSHION  (  SLl'G/SEC  ) 

C   RW"  WATE*  DEnSjTY  (LH*SEC**?/FT»'»'f  ) 

C  Rws=  AIh  DENSlTy  < LB* SEC**2/F T* * <♦  > 

C   TN=  WAVE  HElfiHT  (ft) 

C    U=  FORWARD  <:Pf.-EU  (FT/SFC) 

C    W«  CRAFT  WEIGmT  (LB) 

C   ^F=  FRfQ  OP  wAvES  (RAD/StC) 

C   k'L*  wave  LENGTh  (FT) 

C   I*M"  CRafT  MASS(LB-SEC2/FT  J 

C   WW*  ENCOUNTER  fREQ  FOR  r~EAD  SEAS  (RAD/SFC) 

C   XB=  CUchION  LEnjGTh  (FT  | 

C   Y6=  CUchICN  wIpTh  (FT) 

C  YSw-  SInEWALL  WiDTh  (FT) 

C   Z6=  SEal  hEIght  (FT) 

C  ZIN=  AVr,  WAVE  HEIGHT  UNDER  CUSHION  (FT) 

C  ZOB=  wave  HEir,MT  at  ESOW  (FT) 

C  ZOS=  WAvE  HEIGKj  at  STERN  (FT) 

C   Y<1  )."Zi«POSITIqN  OF  craft 

C   T(2)«F(1)»  VERTICAL  VELOCITv  OF  CRAFT  (FT/SEC) 

C  Y(3)=  air  MASS  IN  CUSHION   (SLUG) 

c  y<5)=pCi,cushiCm  pressure  (LP/<rT2) 

C   F(2)-  ACCELERATION  OF  CRAFT  <FT/SEC**2) 

C   F(3)=  AiR  MasS  FLOW  RATE  (S|.U3/SEC) 
C   Y(9)«  i,\|PUT  wAvE  HEIGHT  (FT) 

C   Y(lif)B  GAP  hFIgHT  FOR  B0w  A\D  STERN  (AVF)  (FT) 
SUBpo'JTInf  f^sim 

COMMON  T,nTJY(2a)iF(?0),STlME*FTIME,HEWDT,lF!  Rt,N 
IF  (nFWOT)  !/3,2 

1  REAn<  3>iC*)     Y»'YBiYSw.ZSW,W>FR/QA,Hl,MPjA^BjTr>ivL»UjPOjAEiPHr 

10  FORMAT  (  «F  1 1*  ,3  ) 

1*  F0RMAT(lFie,2) 

RW»ft?. h 

TN*TN/( 12.0*2*0) 

wM«W/3i?.? 

G«3?.? 

ab»xb*yq 
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TK«53.3<**<*a0»0  +  FR> 

RWS.pA/TK 

WMA«fcM+A 

WF«( p. 0*3,1 ^S*G/WL)**0i5 

W"*  =  wF+(wF*»P«fc')*U/G 

CE*pi,6 

PHAse«3«14?*XB/wL 

R«T\'«YP*wl»pI^(Pi-iASE)/(3.14?*AB) 

F<2  >«^.O000j5i5 

Fl*0.iZ>R0?l/w^A 

F 1 »  71 1 0 

P*«l  »2 

pT»no*ww/(pK*pO)  +  l»pi 

PK3=cO*<?.f*R*S/G)**?»5*AE*0»5/(PO-PA)**0.5 

pK4=*c(",*(?./,*R*S/G)**p.b*2«e*YB*(PO-PA)**0,5 

Cll-RW/PT 

C13»HZ*Rw/PmC 

D12»-PW/PT 

A21=-(R^*XB*YSw*?.0+ftB*Cll)/WhA 

A22»-B/wma 

A23b-AB*Ci3/WMA 

B2l«l»0/WMA 

B22B(RN*XB*YSrt*?,t,"A(^*Dlr)/WMA 

A31«-PK3*clj"PK4*(-l»0  +  Cl  1/RW) 

A33,-<-t<3*ci;*-pK<»*Cl3/R^ 

B32s-P(<L"i».01p-P'<i»»(lt^-«'Diy/Ww() 

WWT»W«*T 

ZOB«TN*SlN(w*T+PHASE) 

ZGS=TN*SIN(WWT-PHASE) 

ZlN«R*SlN(WwT) 

Y(5>«C11*Y(1 )+C13*Y(-j)*Dl2*ZIN 

Y«9)»-TN«sIni(WwT)/0.1 

Y< 1m  )=-Y ( 1 )/0t 10 

Y(  lj  >»-Y(?)/0.5fl 

Y(i3)»«F{?)/10,0 

Y<l3)»Y(3)/(v5 

Y ( 14  > ■ AE/ ( 2 , 0* YB ) * 10 • 0 

Y<l5J»Y<5)/i0»0 

AEB«yB*(Hi-»-(PO   -paj/RW  +  20B-Y(1)  ) 

aES«yB*(H1  +  (P0   -PA)/RW  +  70S-Y(l  )  ) 

IF( AE8.LE.U.0J  AEB=fl,0 

lF<AES.LE,k',ia)  A  E  S  » ;• .  «3 

AE"AgB+AES 

P<3  =  C-n*(?.^*R^S/G)**n»5*AE*t3»5/(PO-PA)**0.5 

A3l*»PK3*Cll-pK**(-l«0*Cl 1 /R W ) 

A33-.PK3*Cl3»PK4*Ci3/RW 

B32s-;;<3*r>l'.>-p<4»'{ltP  +  DlP/Hw) 

F(l >.Y<2) 

F(2)«A21»v(i)+A??»Y(?)+A?3*Y(3)+B21*Fl*B22*ZlN 

F(3)aA31*Y(l )*A33»Y(3»+B?2*ZIN 
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